Systematics of quadrupolar correlation energies 
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We calculate correlation energies associated with the quadrupolar shape degrees of freedom with 
a view to improving the self-consistent mean-field theory of nuclear binding energies. The Generator 
Coordinate Method is employed using mean-field wave functions and the Skyrme SLy4 interaction. 
Systematic results are presented for 605 even-even nuclei of known binding energies, going from mass 
A=16 up to the heaviest known. The correlation energies range from 0.5 to 6.0 MeV in magnitude 
and are rather smooth except for large variations at magic numbers and in light nuclei. Inclusion 
of these correlation energies in the calculated binding energy is found to improve two deficiencies of 
the Skyrme mean field theory. The pure mean field theory has an exaggerated shell effect at neutron 
magic numbers and addition of the correlation energies reduce it. The correlations also explain the 
phenomenon of mutually enhanced magicity, an interaction between neutron and proton shell effects 
that is not explicable in mean field theory. 
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More than 3000 atomic nuclei with different combina- 
tions of proton number Z and neutron number N have 
been experimentally identified so far, and for about 2000 
of them their masses, or, equivalently, their binding en- 
ergies, have been measured [lj. This number is only a 
fraction (about 25%) of all the nuclei which are predicted 
to be bound, i.e. stable against nucleon emission. Since a 
large fraction of the unknown nuclei seems out of reach to 
be studied experimentally, an accurate theory of binding 
energies is needed. For example, one of the important 
motivations for studying nuclei far from stability is to 
better understand nucleosynthesis and the astrophysical 
environment in which it takes place 0. In particular, the 
r-process of heavy-element production leaves an imprint 
of the binding energies of far-unstable nuclei on the mass 
abundances present today. 

Among the theories of nuclear binding energy, the self- 
consistent mean-field theory [3j stands out as a funda- 
mentally justified approach that is also computationally 
tractable over the entire mass table. Global calculations 
of binding energies are now available including both pair- 
ing and deformation effects in the mean field, using en- 
ergy functionals of the nonrelativistic Skyrme form 4] or 
a relativistic mean- field model [5j . A systematic problem 
with such theories is an exaggeration of the increased 
binding at magic numbers, the shell effect. The accu- 
racy of these pure mean field theories is of 1.5-2.0 MeV 
rms errors in binding energies 6] , much poorer than that 
of more phenomenological approaches such as the Finite 
Range Liquid Droplet Model [Tj. 

An obvious route to improve the theory is to include 
correlation effects. Indeed, this has already been done in 
a phenomenological way with good results 0, . Some 
correlation effects can be subsumed in the parametriza- 
tion of the mean-field energy functional, but not all. 
Most obviously, the rotational correlation energy of de- 



formed nuclei is beyond mean-field theory, which can only 
describe the intrinsic state of a rotational band. The re- 
quirements on a theory of correlation energies are severe. 
Besides being computationally tractable, it must be sys- 
tematic, applicable to both spherical and deformed nu- 
clei. One approach, the RPA theory of correlation en- 
ergies, has been applied recently to the binding energy 
problem by Baroni et at, [lOj. However, the RPA is not 
well convergent for short-range interactions, limiting its 
usefulness 

In this letter we report on global calculations of corre- 
lation energies in another well-defined theoretical exten- 
sion of the mean-field approximation |l2j . based on the 
Generator Coordinate Method (GCM). The quadrupolar 
fluctuations about the mean-field solution are determined 
variationally by mixing configurations around the mean- 
field ground state. Our method includes also a projection 
on good angular momentum and particle number and ap- 
plies to all nuclei irrespective of the shape characteristics 
of their ground states. It is thus systematic, which is an 
important criterion for constructing mass tables. It also 
satisfies the goal of including in the correlation energy, 
the rotational energy of deformed states. 

The technical details of our calculation are as follows. 
The self-consistent mean-field equations are solved using 
the method presented in ref. |l3|. Here the wave func- 
tions are represented on a triaxial 3-dimensional spatial 
grid, rather than in an oscillator basis as was done in ref. 
[j, The effective nucleon- nucleon interaction is an 
energy density functional of the Skyrme form; we have 
used the SLy4 parameter set ^ij. This set was fitted to 
nuclear matter properties, nuclear charge radii, as well 
as binding energies of doubly magic nuclei. We also in- 
cluded a local pairing interaction of the form used in ref. 
flij |. but with a reduced strength of vo = —1000 MeV 
fm 3 to compensate the explicit correlations added by the 
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generator coordinate treatment 

The GCM has as ingredients a set of wave func- 
tions generated by constrained mean-field calculations 
and a Hamiltonian for calculating matrix elements among 
those configurations. The constraint that we use is the 
quadrupole operator Q(r) and we minimize the function: 



E[<f>] — X J d 3 rp{r)Q(r). 



(1) 



Here p is the mass density and A is a Lagrange multiplier 
adjusted to get a specific value of the quadrupole moment 
q = J d 3 rp(r)Q(r). The value of q is positive for prolate 
and negative for oblate deformations. The choice of the 
configurations \q) is a computational issue that has to be 
carefully addressed for systematic calculations. If states 
are too closely spaced there is a large redundancy as well 
as an excessive computational cost. We have tailored 
our selection of configurations for a target accuracy of 
0.2 MeV in the correlation energies. In ref. ^(|, it was 
found that this accuracy can be achieved when overlap 
probabilities of neighboring configurations is about 0.5 or 
higher, and we adopt that criterion here. We also restrict 
the deformations to shapes that have energies within 5 
MeV of the mean-field ground state. Typically this gives 
n q = 7 — 20 different configurations in our computational 
space 1171 . 

There is an ambiguity in the choice of the Hamiltonian 
H for use with a mean- field energy functional [lg. We 
use the so-called mixed density for the calculation of the 
density-dependent term in the Skyrme Hamiltonian, as 
is mostly done in the literature |l9j • The GCM energy is 
the minimum of the expression 

(*l#1*> 



where Y$>) is an arbitrary state in the basis of the \q). The 
minimization is performed by solving the corresponding 
discretized Hill- Wheeler equations [20j. The orientation 
of the intrinsic deformed states \q) could be considered a 
continuous degree of freedom to be included in the min- 
imization, but in practice we achieve the same result by 
projecting the intrinsic states on angular momentum zero 
when calculating matrix elements. The theoretical cor- 
relation energy is then obtained in two steps: we deter- 
mine first the energy gain of an angular-momentum pro- 
jected configuration \q) with respect to the mean- field 
minimum, Ej—q, and after that, the additional energy 
gained by mixing configurations with different q magni- 
tudes, Ehw- 

The needed matrix elements (q\q') and (q\H\q') are 
calculated with the computer code described in ref. [l^ . 
Besides performing the rotations to calculate the matrix 
elements needed for angular momentum projection, the 
code carries out a projection of the BCS wave function 
on good particle number N and Z. The number pro- 
jection and the rotation operation are computationally 



very demanding, and it was necessary to find additional 
computational shortcuts to make the global study fea- 
sible. One important savings was the use of a topo- 
logical Gaussian overlap approximation which allows a 
two-point evaluation of the angular momentum projec- 
tion integrals [H IHl. The top-GOA was found to be 
quite adequate, except for very light nuclei, where a 
three-point approximation was sometimes needed. An- 
other place where we made considerable computational 
savings was in the calculation of off-diagonal matrix el- 
ements in the \q) space. In principle, an n-dimensional 
basis requires the calculation of n(n + l)/2 overlaps and 
Hamiltonian matrix elements. These must be done ex- 
plicitly for neighboring configurations, but matrix ele- 
ments for more distant configurations can be estimated 
using another Gaussian Ansatz based on a measure of the 
separation between the configurations I6J. A tricky point 
arises in the mixing of weak prolate and oblate deforma- 
tions which can have very high overlaps. We deal with 
this problem by defining a two-dimensional metric for 
calculating the separation of configurations as described 
in ref. |17| . With these approximations, the computa- 
tional effort for the GCM minimization was reduced by I 
to 2 orders of magnitude, but still requiring on the order 
of 10 17 floating point operations for the entire table of 
605 nuclei. 
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FIG. 1: Typical potential energy landscapes in the deforma- 
tion coordinate. Shown are the curves for 208 Pb, 180 Hg, 170 Hf, 
representing a magic nucleus, a soft nucleus, and a deformed 
nucleus, respectively. See text for explanation. 



We now present three calculations, typical of three dif- 
ferent topologies. In Fig^ are plotted the deformation 
energy curves of a doubly magic nucleus, 208 Pb, a stati- 
cally deformed nucleus 170 Hf, and a nucleus with oblate 
and prolate minima close in energy, 180 Hg. The config- 
uration spaces used are given in Table 1. The 208 Pb is 
a very stiff nucleus, and only needed 4 configurations to 
achieve the target accuracy. In contrast, 180 Hg is a rather 
soft nucleus, and it required 15 configurations ranging 
from -20 b to +32 b to map out the accessible defor- 
mations. The mean-field potential energy landscapes for 
these three nuclei are shown as the dashed lines in Fig.^ 
One sees the narrow minimum at zero deformation for 
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TABLE I: Configuration spaces for typical heavy nuclei. Mass 
quadrupole moments q are given in units of barns. Oblate and 
prolate configurations are listed on separate lines. 



Nucleus 


n q 


q values 


208 Pb 


4 


-10 -5 
+5 +10 


180 Hg 


15 


-20 -16 -14 -10 -6 -4 
+4 +6 +8 +12 +16 +20 +24 +28 +32 


170 Hf 


12 


-20 -16 -13.75 -10 -5 
5 +10 +15 +19.25 +22 +25 +30 



TABLE II: Quadrupolar deformation and correlation energies 
of typical heavy nuclei. Energies are given in MeV. 



Nucleus 


E de f 


Ej = o 


EhW E corr 


208 Pb 


0.0 


-1.7 


0.0 -1.7 


18 °Hg 


-3.0 


-2.5 


-0.5 -3.1 


170 Hf 


-12.2 


-2.9 


-0.4 -3.4 



the magic 208 Pb, a much flatter potential energy surface 
for 180 Hg, and a surface with a strongly deformed mini- 
mum for 170 Hf. The energy gain from angular momentum 
projection is shown by the solid lines in the table. Obvi- 
ously, there is no change for the spherical configurations 
but for finite \q) the projection lowers the energy below 
the mean-field value for every q. Even for 208 Pb, there 
is a net energy gain Ej=q of the lowest projected en- 
ergy with respect to the mean-field ground state. These 
numbers are given in Table 2. 

The second step in our calculation is to apply the GCM 
to mix configurations with different values of q. The ad- 
ditional energy Ehw is given in Table 2 and shown as a 
dot in Fig. 1. The horizontal position of the dot indi- 
cates the average deformation of the mixed wave function 
(see 0| for its precise definition). 

Our calculations include all but the lightest even-even 
nuclei for which the mass is known. The results presented 
in Fig. 2, are displayed as a function of neutron number 
N with isotopes connected by lines. On the top panel, we 
show the energy difference between the SLy4 mean-field 
theory and experiment Here the shell effect is very 
prominent at neutron numbers N = 50, 82, and 126. The 
middle panel shows the calculated quadrupolar correla- 
tions energies. One sees that, for most nuclei, it is be- 
tween 3 and 4 MeV, irrespective of whether the nucleus 
is spherical or deformed. However, near doubly magic 
nuclei, the correlation energy is much smaller. This dif- 
ference has the right sign to mitigate the too-strong shell 
effect of the mean-field theory. The bottom panel in Fig. 
2 shows the binding energies with and without inclusion 
of the correlation energy. There is a distinct improve- 
ment, but one sees that there remains some residual shell 
effect. Since we have only included the axial quadrupo- 
lar field, there could be significant correlation energies 
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FIG. 2: Top panel: Difference of experimental and mean-field 
binding energies for the SLy4 interaction as a function of 
neutron number with lines connecting isotopic chains. Mid- 
dle panel: Theoretical quadrupolar correlation energies of the 
even-even nuclei. Note the expanded energy scale. Bottom 
panel: Difference of experimental and theoretical binding en- 
ergies including the correlation energy. Comparison between 
binding energies obtained with and without the correlation 
contribution. 

arising from other kinds of deformation or from pairing 
vibrations around closed shells [lfj. Also, it is clear that 
a more complete theory requires a refitting of the energy 
functional to take into account the contribution from the 
correlation energy. 
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FIG. 3: Two-proton gaps for Pb and Sn isotopic chains. 

A feature of the binding systematics that has been dif- 
ficult to understand in mean-field theory is "mutually 
enhanced magicity" |23( , the enhancement of shell effects 



when both neutrons and protons form closed shells. This 
behavior can be seen in the neutron dependences of the 
two-proton gap, defined as 

S 2p (N,Z) =E(N,Z-2)-2E(N,Z)+E(N,Z+2). (2) 

In the independent-particle shell model, &2 P is propor- 
tional to the difference between the Fermi energies of 
nuclei differing by two protons and thus it is a measure 
of shell gaps. This quantity is plotted in Fig. 3 for the 
Z = 50 (Sn) and Z = 82 (Pb) proton gaps, with the 
experimental points shown by black diamonds. There 
is clearly a reduction of bi v when going away from the 
doubly magic 132 Sn with 82 neutrons and doubly magic 
208 Pb with 126 neutrons. Neglecting deformations, the 
mean- field prediction for bi v is quite flat, as shown by the 
short-dashed lines on the graphs. This reflects the inde- 
pendence of the single-particle proton gap on the neu- 
trons in that approximation)^. Allowing deformations 
in the mean-field energies already gives a change into the 
right direction, as shown by the long-dashed lines. The 
calculation including the full quadrupole correlation en- 
ergy is shown by the solid lines. One sees that it brings 
the calculations even closer to experiment. The corre- 
lation energy thus provides a plausible explanation of 
mutually enhanced magicity. This conclusion is corrobo- 
rated by a recent calculation of bi v in the Sn chain using 
the Bohr Hamiltonian, showing also the large improve- 
ment brought by quadrupolar correlations |24j . 

In conclusion, we find that quadrupole correlations be- 
yond the mean field capture a significant part of the bind- 
ing energy that is not describable in mean-field theory. 
We have shown that global calculations of these corre- 
lations are today feasible and could be incorporated in 
the future in the construction of effective interactions. 
Further improvements should be obtained by including 
other correlation effects related to octupole deformations 
or pairing vibration 0, . Our basis does not contain 
triaxial shapes, but we do not believe they play as impor- 
tant a role. One of the major challenges which remains 
open is to include neutron-proton pairing effects, which 
may be responsible for the scatter of the residuals around 
the light N = Z nuclei in Fig. 2. 
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